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We propose the reahzation of topological quantum states with cold atoms trapped in an optical lattice. We dis- 
cuss an experimental setup that generates a two-dimensional hexagonal lattice in the presence of a light-induced 
periodic vector potential, which represents a realization of the Haldane model with cold atoms. We determine 
theoretically the conditions necessary for observing the topological states and show that two of the key condi- 
tions are: 1) the realization of sharp boundaries and 2) the minimization of any smoothly varying component of 
the confining potential. We argue that, unlike their condensed matter counterparts, cold atom topological quan- 
tum states can be i) "seen", by mapping out the characteristic chiral edge states, and ii) controlled, by controlling 
the periodic vector potential and the properties of the confining potential. 



A solid state insulator can be defined as a system with local 
electronic properties 11]. Consequently, insulators are insensi- 
tive to boundary conditions. Standard band insulators, which 
are characterized by the existence of a bulk energy gap, satisfy 
this definition and are topologically equivalent, as they can be 
adiabatically transformed into each other without crossing a 
phase transition. However, the existence of a bulk gap is not a 
sufi&cient condition to insure the locality of all electronic prop- 
erties and the insensitivity to boundary conditions. Certain 
strongly correlated systems, such as the fractional quantum 
Hall fluids, ofifer examples of phases having bulk gaps but be- 
ing topologically distinct O. Such topological insulators can 
exist even in the absence of interactions and typical examples 
are the integer quantum Hall fluids or the quantum spin Hall 
states El ^4, 5 1 and their three dimensional generalizations 
One defining characteristic of these systems is the existence of 
chiral gapless edge states that are robust against disorder ef- 
fects and interactions. The existence of these edge states and 
their basic features are intrinsically linked to the topological 
properties of the system. However, their detailed structure is 
dictated by the boundary. Controlling the boundary properties 
is a rather difl&cult task for solid state systems, but it could be 
realized for ultracold atoms in optical lattices. 

In this work we propose the realization of a topological in- 
sulator in an optical lattice. More importantly, we establish 
that optical lattices allow for the existence of a "topologi- 
cal metal" characterized by a chiral boundary-induced edge 
mode, although the bulk is metallic. A great advantage of the 
optical lattices is that various terms in the Hamiltonian as well 
as the boundary conditions can be explicitly controlled exper- 
imentally, which in contrast to solid state systems, allows for 
tuning the properties of the edge states. The realization of a 
topological insulator with cold atoms opens a series of very 
exciting prospects: i) the possibility of a direct observation of 
the edge states |7 |, the hallmark of the topological insulating 
phase, ii) the possibility of testing the stability of the chiral 
edge modes in the presence weak disorder and interactions, 
and iii) the possibility of studying transitions between a topo- 
logical insulator and other phases. On the other hand, the main 
challenges in building a topological insulator with cold atoms 
are: i) generating the vector potential, ii) controlling the trap 



potential and manipulating the boundaries, and iii) measuring 
a topological insulating state, i.e., imaging the edge states. 

To study the boundary eflTects in a topological quantum 
state, we propose the realization of the Haldane model HI 
in an optical lattice. The trapping potential is given by the 
superposition of three co-planar standing waves character- 
ized by the wave- vectors ki = (0, |^), k2 = (-^. ^) and 
ks = (--^, ^), respectively. The minima of this poten- 
tial generate a hexagonal lattice with lattice constant a. The 
crucial ingredient of the setup is an eflfective vector poten- 
tial A(r) that generates a periodic "magnetic" field with zero 
total flux trough the unit cell. We propose the use of a 
light-induced gauge potential that can be realized in a sys- 
tem of multi-level atoms interacting with a spatially modu- 
lated laser field |[8l|9l[l0l[ni[I2[l3l[Il|T5l|Ti[I3 Within 
these schemes, the multi-level atoms interact with laser beams 
characterized by spatially varying Rabi frequencies and ex- 
perience an eflTective pseudo-spin dependent gauge potential. 
Since our proposal does not require spin-dependent gauge po- 
tentials, one may even be able to utilize a simpler scheme ifTSl . 
The eflfective single particle Hamiltonian is 

^ 3 
= ^ (P - A(r))2 + Vo 2 cos2(k,r) + y,(r), (1) 

i-l 

where m is the particle mass, p the momentum and A the 
vector potential. The second term in Eq. ([T]) generates 
the two-dimensional hexagonal optical lattice with lattice 
constant a, while the last term, V^, contains additional con- 
fining terms that determine the properties of the boundaries. 
The role of will be discussed in detail below. A simple 
vector potential that generates zero "magnetic" flux trough 
a unit ceU is A(r) = aJi(r) with Ji^(r) = sin(^^) and 
Jiy = 0, where a is the strength of the gauge potential. 
Notice however that the term A? 1 2m from Eq. ^ represents 
an extra contribution to the lattice potential that does not 
have hexagonal symmetry and, therefore, will distort the 
lattice. This distortion does not aflfect the nature and the basic 
properties of the topological insulator if a does not exceed 
a certain critical value, but has quantitative implications for 
the band structure. To avoid the formal complications of 
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FIG. 1 : Optical lattice potential formed by the superposition of three 
standing waves and generating a two-dimensional hexagonal lattice 
with lattice constant a. The effective confining potential along the 
segment (0, 0) - (0, 2a) (red line in the upper right panel) is shown 
in the absence of a vector potential (blue circles) and for A 9^ (red 
triangles). Inset: typical cluster used in the calculations consisting of 
a disc-shaped piece of hexagonal lattice with radius R ^ 39a. 

dealing with this Of-dependent lattice distortion, we use in 
our calculations a symmetrized vector potential with J{(y) = 
[sm(4ny/3) + cos(27rx/ VS) sin(27r///3), - Vs sin(27rjc/ VS) 
X cos(27rz//3)]. The fact that A? 1 2m acts as an effective lat- 
tice potential suggests that ideally the trapping and vector 
potentials could be realized using the same set of laser fields. 

We solve the quantum problem associated with Hamil- 
tonian ([T]) within a simple tight binding approximation. 
Throughout the paper we use the recoil energy Er = ^^^^ 
as energy unit and the lattice constant a as length unit. 
The total effective lattice potential Viatt^ which includes the 
term A? 12m, with minima at the nodes of a hexagonal lat- 
tice (see Fig. [T]) has near these minima the form Viau ~ 
majl/2(Sx^ + Sy^). This suggests the use of the s-wave or- 
bitals 0Q^(r) = ^/2/(7lc) Qxp[-(r - r/)^/c] as a possible simple 
basis for the tight-binding approximation. Here represents 
the position of a lattice site, c = (4Er)/(n:^oJo)a^ and we have 
Vo = l2Er/7r'^(a^/c^-7i^a^a^/4). The approximation holds as 
long as the s-band is well separated from the p-bands, which 
is the case for c < 0.25, i.e., for deep enough optical lattices. 
However, because the the second-neighbor hopping is crucial 
for generating the topological states 14], a small value of c 
will make this effect practically unobservable. In this study 
we choose c = 0.2. The other independent parameter is a, 
which in the calculations will be either a = (zero vector 
potential) or a = 2. The hopping parameters for the effective 
tight-binding model are given by tij = The key 

contributions coming from the vector potential, (0Q^|pA|0Q ^), 
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FIG. 2: Finite size equivalent of the Hofstadter "moth" for a cluster 
with a hard-wall boundary in the presence of a periodic vector po- 
tential with Of = 2. The coordinates of each dot represent the orbital 
momentum (x-axis) and the energy (y-axis) of a particular state. The 
"edge" vs. "bulk" character of each state is shown by the color code, 
which represents the relative boundary contribution to the norm, y„ 
(see main text), ranges from 1 for edge states (blue) to for purely 
bulk states (red). Notice the chiral nature of the edge states mode that 
populates the bulk gap of the insulator. The corresponding density of 
states (DOS) in shown in panel (2). For comparison, in panel (3) we 
show the DOS in the absence of a gauge potential. 

vanish if i and j are nearest neighbors and are non-zero for 
second-order neighbors, generating a chiral contribution to the 
Hamiltonian. Finally, because the s-wave orbitals are not or- 
thogonal, we have to determine their overlap matrix elements 
The resulting tight-binding problem is solved nu- 
merically for clusters containing up to 3696 sites (see Fig[T]). 

We start by solving the tight-binding problem for a finite- 
size cluster consisting of a disc with hard-wall boundary con- 
ditions, Vc(r) = 00 if r > R and Vdr) = if r < R. The density 
of states (DOS) for this system is shown in Fig. |2] (panels 2 
and 3). In the absence of a vector potential (a = - panel 
3), this quantity is similar to the graphene DOS. The charac- 
teristic V-shaped gap associated in the infinite lattice with the 
existence of the Dirac points is weakly distorted at low ener- 
gies by finite size effects that become negligible as the cluster 
size increases. In the presence of a vector potential with a = 2 
a finite size gap opens in the density of states (panel (2) Fig. 
|2]). However the DOS in the "gap" is not exactly zero. To de- 
termine the nature of the residual in-gap states, we calculate 
the orbital momentum of each single-particle state </^„(r) and 
the relative contribution to the norm {il/nWn) coming from 
a narrow ring ?>la < r < 39a positioned at the edge of the 
system. This contribution vanishes for bulk-type states and is 
of order one for edge states, i.e., it represents a measure of the 
edge-type character of a given state. The corresponding spec- 
trum is shown in Fig. [2] panel (1). This picture is the finite 
size equivalent of the Hofstadter "moth" for a hexagonal lat- 
tice and a periodic "magnetic" field with zero flux through the 




FIG. 3: Left panels: Contour plots of p„(r) = |<An(r)p for the states 
marked in Fig. |3] The quantity p„(r) is the product between a com- 
mon factor associated with the underlying hexagonal lattice struc- 
ture and a state-dependent envelope function. The corresponding en- 
velope functions are shown in the right panels. (A) represents the 
ground state, (B) is a typical edge state and (C) is the lowest energy 
bulk state from the upper band. The weak edge contributions in (C) 
are due to finite size eff'ects and vanish in the large cluster limit. 



unit cell. The coordinates of each dot represent the orbital mo- 
mentum (x-axis) and the energy (y-axis) of a particular state. 
The edge vs. bulk character of the state, as quantified by the 
parameter is revealed by the color code: blue for edge states 
and red for bulk states. As evident from Fig. |2] the spectrum 
is characterized by a gap for the bulk (red) states. Within this 
gap, there is an chiral edge mode (blue states). The chirality 
of the edge mode, i.e., the sign of its orbital momentum, is 
determined by the sign of a. Note that the sum of the orbital 
momenta of all the single-particle states is identically zero for 
any value of a. To have a spatial characterization of the single- 
particle quantum states, we show in Fig. [3] the contour plots 
of p„(r) = |<A«(r)p for several states marked in Fig. [2j A - 
the ground state, B - a typical edge state and C - the lowest 
energy bulk state from the upper band. Note that each density 
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FIG. 4: System with soft boundaries. The confining potential is 
shown in the top panels for: a) quartic wall, b) linear step, and c) 
harmonic trap plus infinite wall. The turquoise region corresponds 
to energies smaller than the gap Aq,, while in the yellow region the 
energy is smaller than the bandwidth W. A confining potential with 
characteristic length much smaller than the size of the system,(a) and 
(b), preserves the bulk gap (topological insulator). The in-gap fea- 
tures appearing in the density of states (middle panels) are all due to 
edge states (see also Fig. |5]). In a smoothly varying confining po- 
tential, case (c), the bulk gap collapses (topological metal, see lower 
panel). In cases (a) and (b) the lower band is unaff'ected by the details 
of the confining potential, while the in-gap structures (middle panels) 
are similar. The rapid oscillations at low-energies in case (c) indicate 
the formation of harmonic oscillator levels. 

function p„(r) consists of a common factor 0o^r ~ ^i)] as- 
sociated with the underlying lattice structure multiplied by a 
specific envelope function. The structure of the envelope func- 
tions is shown in the left panels of Fig. [3] 

Next, we address the key question concerning the role of 
the confining potential Vc and the dependence of the spec- 
trum on the boundary conditions. We replace the hard wall 
boundary by a) a quartic wall, b) a linear step potential, and 
c) a harmonic potential plus a hard wall (see Fig. |4] top 
panel). The corresponding expressions of the confining po- 
tential are V^(r) = A^(r - Rq)^ if r > Rq (and otherwise), 
V^(r) = Min [4(r - Ri)/(R2 - Ri), 4] if ^ > (and oth- 
erwise), and V^(r) = A^r^ if r < R2 (and 00 otherwise), re- 
spectively. To define a characteristic length- scale associated 
with the confining potential, we introduce the radii R^, such 
that Vc(Ra) = Aa and Rw, with the property VdRw) = ^• 
Here is the bulk gap and W the sum of the bandwidths 
of the lower and upper bands for a system with hard walls. 
The relevant length scale for a soft boundary produced by the 
confining potential Vc is given by dc = Rw - Ra- Our nu- 
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FIG. 5: In-gap states for a system with soft boundary (quartic wall). 
The color code is the same as in figure Fig [2] Note that all the in-gap 
states have edge -like character. In addition to the typical edge mode 
consisting of chiral edge states (A), notice the presence of Tamm-like 
edge states (B) with opposite orbital momentum. 

merical calculations show that a topological insulator can be 
realized provided dc ^ Ra, i.e., the boundary has a charac- 
teristic length much smaller than the size of the system. For 
example, cases (a) and (b) in Fig. [5]correspond to ^ 32.5a 
and dc ~ 3.5a and in both cases the gap for bulk states is pre- 
served. However, in contrast with hard wall case characterized 
by a featureless residual in-gap DOS (see Fig. |2] panel 2), a 
system with soft boundaries has a nontrivial structure of the 
residual DOS (Fig. [4j middle panels). This structure emerges 
from two causes: 1) the orbital momentum of the chiral edge 
mode acquires a more complicated energy dependence, and 
2) additional edge states, that do not belong to the chiral edge 
mode, develop inside the gap. Both these points are illustrated 
in Fig. [5] The B-type edge states are Tamm-like states, which 
are formed due to the rapid variation of the confining potential 
and are not related to the topological properties of the insula- 
tor. Finally, in the presence of a confining potential with a 
smoothly varying component (case (c) in Fig. |4]) the bulk gap 
collapses and the system becomes metallic. Note that, in the 
presence of the periodic vector potential A(r), the system still 
has chiral edge states even in the metallic phase, a situation 
similar to the existence of surface states in doped semiconduc- 
tors in the presence of spin-orbit coupling 1 19 |. We conclude 
that, in order to realize a topological insulator with cold atoms, 
one needs to produce a sharp boundary and to minimize any 
smoothly varying component of the confining potential. 

Probing a topological state in condensed matter systems in- 
volves transport measurements. This would be a rather dif- 
ficult task for atomic systems in optical lattices. Since the 



nontrivial topological properties of a system are a feature of 
the single particle Hamiltonian which is best revealed by the 
presence of the chiral edge states, we propose here the direct 
observation of these edge states in cold atomic systems, some- 
thing one cannot easily realize in the condensed matter con- 
text. This involves two steps: 1) loading bosons into the edge 
states and 2) imaging the atoms. Initially, the optical lattice 
is loaded with atoms and cooled, so that the bosons occupy 
only the lowest energy single particle states. Next, angular 
momentum is transferred to the atoms using, for example, a 
2-photon stimulated Raman process |20|, thereby populating 
states of higher energy and angular momentum comparable to 
that of the edge states. Finally, one can use resonant tt pulses 
to further excite the atoms into the edge states. To image the 
atoms, one can use a direct in-situ imaging technique 1211 . 
Although this technique is able, in principle, to discriminate 
between individual atoms in an optical lattice, a much lower 
resolution is required for imaging the edge states. 

We propose the realization of topological quantum states 
with cold atoms in an optical lattice. A combination of a 
hexagonal optical lattice potential and a periodic light-induced 
vector potential represents a cold-atom realization of the Hal- 
dane model. We show that such a system is characterized by 
chiral edge modes, which are the signature of a topological 
quantum state. Observing these edge sates is an eflTective way 
of "seeing" a topological phase. We find that necessary condi- 
tions for realizing a topological insulator with optical lattices 
are the realization of a sharp boundary and the minimization 
of any smoothly varying component of the confining poten- 
tial, e.g., of the harmonic confining potential. Controlling the 
confining potential opens the possibility of testing the stability 
of the chiral edge modes against weak disorder and interac- 
tions and, together with control of the vector potential, offers 
a knob for tuning the system from a topological insulator state 
to a standard insulator or a metallic phase. 
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